Distributed Delays Facilitate Amplitude Death of Coupled Oscillators 
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Coupled oscillators are shown to experience amplitude death for a much larger set of parameter 
values when they are connected with time delays distributed over an interval rather than concen- 
trated at a point. Distributed delays enlarge and merge death islands in the parameter space. 
Furthermore, when the variance of the distribution is larger than a threshold the death region be- 
comes unbounded and amplitude death can occur for any average value of delay. These phenomena 
are observed even with a small spread of delays, for different distribution functions, and an arbitrary 
number of oscillators. (http://link.aps.org/abstract/PRL/v91/e094101l 
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Coupled oscillators constitute an effective and popu- 
lar paradigm for the study of interacting oscillatory pro- 
cesses in the physical and biological sciences jl|. The 
rich dynamics arising from the interaction of simple units 
have been a source of interest for scientists modeling the 
collective behavior of real-life systems. Among the most 
widely studied phenomena is synchronization, where in- 
dividual units oscillate at a common frequency and phase 
when coupled 2] . Synchronization may be observed even 
under weak coupling; so it has usually been studied 
through reduced models that retain only the phase in- 
formation along the limit cycles. With stronger coupling 
further interesting behavior is possible, whose investiga- 
tion requires the use of full models that include the ampli- 
tudes of the oscillators. An example is amplitude death, 
which refers to the quenching of oscillations under cou- 
pling, as the system evolves to a stable equilibrium 
If the information flow through the coupled system is in- 
stantaneous, amplitude death occurs when the individual 
oscillators have sufficiently different frequencies. 5, 6, 7]. 
On the other hand, if the information from one oscilla- 
tor reaches the others after a certain time delay, which 
may be due to finite propagation or information process- 
ing speeds, then even identical oscillators can experience 
amplitude death when coupled Recent experimental 
and theoretical studies have confirmed the role of delays 
in inducing amplitude death 0, 0, 0| . 

While the importance of time delays in amplitude 
death is now clear, studies in this area have so far been 
confined only to discrete, or constant, delays. In other 
words, it has been assumed that information reaches from 
one unit to another after a fixed time t which is unchang- 
ing as the system evolves, and moreover, the units act 
only on the instantaneous value of the received informa- 
tion and forget any previous values. Such discrete-delay 
models often fail to adequately describe physical systems 
by neglecting the possibilities that (a) the quantity r may 
only be approximately known, (b) it may only represent 
an average value of a quantity that varies between pairs 
of oscillators in a network or (c) varies in time through a 
process involving unmodelled factors, and (d) the oscilla- 



Copyright: American Physical Society 

tors may incorporate "memory" effects by using the past 
history of the received information. The first possibility 
is certainly an issue in any experimental situation, (b) is 
typical when considering large arrays, and (c) and (d) are 
particularly significant in biology and neurology. Because 
of these shortcomings, models based on distributed delays 
have been proposed as early as the time of Volterra 
and used in such areas as biology [l3| . ecology 
neurology |Q| , viscoelasticity , and economics 
has especially been pointed out in the biological sciences 
that distributed delays lead to more realistic models 1^ . 
In this Letter we consider the effects of distributed delays 
on amplitude death, and show that even a small spread 
in the delay distribution can greatly enlarge the set of 
parameters for which amplitude death occurs. 
The oscillators studied here are described by 



Z{t) = (1 + icoo- |Z(i)|2)Z(t), 



(1) 



where Z{t) is a complex number and ojq > 0. Equation 
([T|) represents the normal form for a supercritical Hopf 
bifurcation, and has been used to describe limit cycle 
oscillators where oscillations arise through such a bifur- 
cation. It has an unstable equilibrium at zero, and an 
attracting limit cycle Z{t) — exp(ia;oO with frequency 
lhq. Suppose a pair of such oscillators are coupled with 
general time delays: 



Zi(i) = (1+ic.i- izi(t)nzi(t) 



K 



f{r')Z2{t-T')dr' -Z,{t) 



Z2{t)^il+iLU2-\Z2{tf)Z2{t) 



f{T')Zi{t~r')dr' -Z2{t) 



(2) 



(3) 



Here, K is a, number quantifying the strength of coupling, 
and / represents a distribution of delay values. When / is 
the delta function 6{0) one obtains the system considered 
in Q, where the oscillators interact without delay. Sim- 
ilarly, the choice 6{t) with t > gives the system with 
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a discrete delay which was studied in [g. In the general 
case / is a probability density over an appropriate inter- 
val, which addresses the shortcomings mentioned in the 
above paragraph. We shall show that the variance of / 
has a significant effect on the dynamics of the system. 

When the system experiences amplitude death 

its zero solution becomes stable. To investigate the 
stability, the system is linearized about zero. The 
characteristic equation is found by making the ansatz 
{Zi{t), Z2{t)) = Cexp(Ai), ^ e and is given by 

{l + iuji^K~\){l + iL02-K-\)-K'^[F{\)f =0 (4) 

where F is the Laplace transform of /. For definiteness, 
the analysis of Q will be illustrated for uniformly dis- 
tributed delays over the interval r±a, i.e. /(r') = l/(2a) 
if |r — r'l < a and zero otherwise. However, the idea is 
the same for other types of distributions. We first con- 
sider oscillators with identical frequencies, wi = a;2 = wq. 
This is the more stringent case for stability, because dif- 
ferent frequencies can stabilize coupled oscillators even in 
the absence of delays whereas identical frequencies can- 
not @ . The analysis is based on the observation that 
as parameters are varied the stability of the origin may 
change only if an eigenvalue A crosses the imaginary axis. 
In this critical situation X = iu for some real w, and by 

m 

{l+i{ujo-^)-Kf -K^l'^e-^''^^ (5) 

where 

, . \ siTi{Loa)/{uja) if Loa^O, 
^ = ^("'") = | 1 if ..a = 0. 

Separating ((SJ into real and imaginary parts and rear- 
ranging yields 

{l-j^)K^-2K + l = -{uj-ujo'f (7) 
tan(c.r) = (8) 

This pair of equations describe a set of parametric curves 
on the T-K plane in the parameter cu. For each value of 
oj, K is found from ||7J|, and substitution into (jHJ gives 
the corresponding values for t. The stability region is 
determined by computing the critical curves Q-JSJl, and 
following the direction of movement of the purely imag- 
inary eigenvalues as parameters are varied. The latter 
information is obtained from the quantities Ile{dX/dK) 
and Ke{dX/dT) calculated from Q by implicit differen- 
tiation on the critical curves. For the parameter values 
obtained by this procedure, amplitude death is indepen- 
dently confirmed by numerical simulation of the coupled 
system Q-©- 

The stability region in the K-t parameter plane is 
shown in Fig. ^ for uniformly distributed delays over 
T ± a, with T > a so that the delays involved are not 
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FIG. 1: Enlargement of the stability region as the spread of 
the delay distribution increases: (a) a = 0, (b) a = 0.007, (c) 
a = 0.008, and (d) a = 0.02. The value of loq is 30. 

negative. As a 0, the uniform distribution approaches 
S(t); hence a serves as a parameter to compare the dis- 
crete and the uniformly distributed delays having the 
same mean value t. This is equivalent to quantifying 
the effects of the distribution by its standard deviation 
a = a/y/3. When the delay is discrete the stability re- 
gion consists of three disjoint and bounded sets (Fig.QjL) 
which deform continuously as a is increased from zero 
(Fig. ^). The enclosed area increases with the vari- 
ance of the distribution, and at some critical value the 
stability region becomes unbounded in the r-direction 
(Fig. ^). The critical value of a for this qualitative 
transition depends on luq and it can be very small; it 
is slightly below 0.008 when luq = 30, corresponding to 
a — 4.62 X 10"^. Increasing a further results in a con- 
nected and unbounded set of parameter values for ampli- 
tude death (Fig.^), which persists for all larger values 
of a. At this stage, there is a large interval of values for 
the coupling strength K which causes amplitude death 
regardless of the mean value t. By contrast, the discrete 
delay can cause amplitude death only for a very limited 
range of delay values (Fig. ^). Note that the ratio cr/r 
of the standard deviation to the mean of the distribution 
can be quite small, showing that even a relatively small 
spread in delays may induce amplitude death. 

Distributed delays further facilitate amplitude death 
through the parameter ujq. The stability region in the 
parameter plane depends on the value of wq, and there 
exists a minimum value ujmin such that if coq < ujmin then 
amplitude death does not occur for any K oi t. A value 
of uJmin = 4.812 has been reported for discrete delays 
IQ. Distributed delays can induce amplitude death even 
when Wo is smaller than this value, as Fig. |21 shows. The 




FIG. 2: Minimum value of loq for which amplitude death is 
possible decreases as the spread of the delay distribution in- 
creases. 
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FIG. 3: Amplitudes of the limit cycles of individual oscillators 
versus the standard deviation of the delay distribution, for 
common distribution functions with the same mean value r — 
0.5. Other parameters are K = 30 and ujq = 30. 



curve has been numerically calculated by decreasing the 
value of ujQ at a fixed cr for the uniform distribution until 
the stability region disappears. 

The local stability change caused by distributed delays 
is reflected in the global dynamics by the annihilation of 
the limit cycles of the oscillators. Fig. 13 shows the am- 
plitudes of the oscillators, which gradually decrease with 
increasing standard deviation of the delay distribution. 
Near a = 0.0085 the amplitude becomes zero as the limit 
cycle collapses to the origin and the amplitude death sets 
in. Numerical simulations with random initial conditions 
indicate that for a beyond this value the origin is the only 
attractor for the coupled system. Notably, this behavior 
is largely independent of the particular distribution cho- 
sen for the delays, as seen in Fig.|2| 

When the oscillators have different intrinsic frequen- 
cies the situation is similar, except that amplitude death 
occurs for even a larger parameter set, since a large fre- 
quency difference by itself is known to cause death. Fig.H 
compares the stability regions for discrete and distributed 
delays. As the spread of the delays increase, the stability 
region is enlarged and extended towards the horizontal 
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FIG. 4: Stability regions for coupled oscillators with different 
frequencies when the delay is (a) discrete at r = 0.5, and (b) 
uniformly distributed over 0.5 ± 0.02. The mean frequency of 
the oscillators is fixed at 30. 



axis, so amplitude death becomes possible also for a small 
(or zero) frequency difference. 

The stabilizing effects of distributed delays carry over 
to any number of oscillators. This is illustrated in the 
following system of N globally coupled oscillators: 



Z,{t)^il + i^,-\Z,it)\^)Z,{t) 



2K^ 

N 2^ 



k=l 



f{T')Zk{t-T') dT'-Z,{t) 



(9) 



where j = 1, . . . , TV. In this model the use of distributed 
delays is further motivated by the fact that in a physi- 
cal network the delays between units cannot be expected 
to be identical but can be more naturally described by 
a probability distribution for large N . For N = 2 

reduces to Q^©- The system was studied in 

with no delays, and in ^ with discrete delays. For iden 
tical oscillators an analysis similar to above yields the 
region of amplitude death for @. The limiting shape of 
this region as A'^ — > oo is depicted in Fig. O for uniformly 
distributed delays and LOj — 10 for all j. As before, dis- 
tributed delays enlarge the stability region, and there 
exist values of K for which amplitude death occurs re- 
gardless of the mean value r of the delays. (The stability 
region is enlarged further if the frequencies ojj are not 
identical, similar to the case shown in Fig. HI) It is in- 
teresting to compare the results to those in ^Sb l3] j where 
it was shown that a sufficiently large spread in the fre- 
quencies LHj can cause amplitude death. Here a similar 
conclusion holds for a spread in the delays. 

A detailed mathematical analysis of stability under 
distributed delays is too lengthy to include here; how- 
ever, a brief description of the basic ideas will help clarify 
the role of delay distributions. Thus consider the char- 
acteristic equation Q. If = the eigenvalues are 
A = 1 — K + iuoj, j — 1,2, so the system is asymptotically 
stable ior K > 1. It follows that stability can also be 
achieved for |F| sufficiently small. Letting X = (3 + iuj, 



FiX) 



f{r')e 



'^'^ [cos lot' — i sin wt'I 



dr' 
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Finally, distributed delays give a natural way to model 
memory effects in interacting systems. They are thus 
particularly suitable in fields of neuroscience, cognition, 
and the general analysis of complex systems. Delays are 
expected to be a source of further interesting results in 
these highly active areas of investigation. 



FIG. 5: The limiting shape of the stability region for an in- 
finite number of globally coupled oscillators, calculated for 
delays uniformly distributed over r±0.07. The area enclosed 
by the dashed line is the stability region for the discrete delay. 



which can be interpreted as a (weighted) average of the 
quantity in brackets over an interval determined by /. If 
w 7^ the integrand is oscillatory about zero, and smaller 
values of |i^| may be obtained if the average is taken over 
a larger interval, i.e. if / has a large variance. It can 
be assured that lo is bounded away from zero provided 
K is not too large. In this way, an interval of values of 
K is obtained for which the origin is stable. A rigor- 
ous argument involves assuming /? > and obtaining a 
contradiction resulting from \F\ being small. It is also 
intuitively plausible that if / has a sufficiently large vari- 
ance, then |F| is small regardless of the mean value or the 
precise shape of /. This explains Figures^ and[51 which 
show stability for a range of coupling strengths and arbi- 
trary mean value of delay, as well as Fig. |31 which shows 
similar behavior for different distribution functions. 

The importance of amplitude death has been noted by 
many authors in relation to various physical and biolog- 
ical phenomena, ranging from Belousov-Zhabotinskii re- 
actions to cardiac arrhythmias 4,0,@|. For instance, the 
cessation of rhythmic activity in biological systems may 
be related to certain pathologies. Our findings suggest 
that in certain cases the variance of the delays, rather 
than their average value, could be the relevant quantity 
responsible for the quenching of oscillations. Delay distri- 
butions has a stabilizing effect on the interconnected sys- 
tem, similar to that of frequency distributions treated in 
previous works 0, . This implies that amplitude death 
is a rather common and robust dynamical behavior for 
interacting oscillatory processes, since real-life networks 
inevitably involve variances in both the frequencies and 
the connection delays. From a practical point, the prop- 
erties of distributed delays are expected to be helpful in 
modelling observed phenomena. For instance, a situa- 
tion where amplitude death is experimentally observed 
but not predicted by an existing model might imply the 
presence of distributed delays in the physical system. In- 
troducing a small variance in delays may in many cases 
lead to a better reconciliation of theory with experiments. 
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